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Abstract 


This report (Preface PI, P2; Parts A-E) is essentially a 
Summary of the five M.S. Theses and Ph.D. Dissertations completed 
during the contract year for contract No. DA-49-146-XZ-494. The 
titles of these theses and dissertations are the heading titles 
for each of the five Parts A-E. 

This Final Report is more than a mere Summary since all results 
have been condensed, compared, corrected and interpreted as a group. 
Figures A8, C14 and D4 exhibit count-rates and/or particle flux 
densities, both for the Natural and the ARGUS Events 1 and 2 radia¬ 
tions, in 3-dimensional B-L coordinates. These 3-dimensional plots 
exhibit most clearly the coordination of all the data and are pre¬ 
sented for the first time in this Summary Final Report. 



PI. Preface - Glossary 


SYMBOL 


DEFINITION 


A Nucleon number (Ref. PI, Sec. 1) 

A 

A Vector, along axis of minimum moment of inertia Ex¬ 

plorer 4. Also the direction of the satellite symmetry 
axis. (Ref. Pl } Fig. 2.2). 

A e ^ Effective area (product of counter efficiency and geo¬ 

metric area A) of directional counter Explorer 4 (cm ) 


a. 

l 


Scale height for atmospheric model defined by 

1 N i 

In r;- , where 


a. = 

l 


h. h. 

l+l l 


N 


i+1 


a 


A 

P 


A 

B 


A 


B 

o 


B 


m 


Bj. 


N. - number/unit volume at altitude h. where h. is the 
altitude in kilometers (Ref. P2) 1 1 

Equatorial radius of international ellipsoid (km) 

Daily equivalent planetary amplitude (Ref. P3, Sec. 1.5) 

Three-hour equivalent planetary amplitude (Ref. P3, 

Sec. 1.5) 

Geomagnetic flux density vector (gauss) 

Geomagnetic flux density unit vector, B/B 

Magnetic flux density unit vector at the time the satellite 
inertial coordinate system is established (time of any rf- 
null prior to data to be analyzed) 

Mirror point value of B (gauss) 

A vector in the plane perpendicular to the magnetic flux 
density vector B (gauss) 


B v ,B w ,B^ Components of a unit vector B in Vernal Equinox coordi- 

1 nates 

B e Equatorial value of B 

b Polar radius cf international ellipsoid 

A 

b.,b 0 ,b Components of B on the X s> Y s ,Z s axes of satellite 

inertial coordinate systlm (Ref. PI, Eq. (3.14)) 


Pl-1 



c 

c 

c 

w e > 

^obs 

C(6) 

c(?) 

E 

a 

E 

e 

f(h') 

F(oL,<f))J(l 

(f) 

H 

H ’ H R 

h 

I 

1(E) 


i(H) 

I M 
P 


Normalized (divided by B) components in topcentric co¬ 
ordinate system of ephemcris (north, east, down) 

Omnidirectional count-rate corrected for deadtime (sec j 

Speed of light in free space 

Observed count-rate corrected for counter dead-time (sec"*) 

Observed directional count-rate for detector A (sec - *) 

Omnidirectional observed count-rate (sec -1 ) 

True directional count-rate (0^(0) c^rregted for deadtime, 
efficiency and effective area) (see - cm ) 

Auto-covariance function (Ref. PI, P4, Eq. (A1.2)) 

Kinetic energy of a particle 

Electric field vector 

Electronic charge 

Mirror-point density (Ref. P2, Eq. (2.24)) 

Fraction of particles with directions within dCL at 
(Ref. P3, Sec. 4.S) 

Fraction of fission electrons that penetrate a shielding 
(Ref. P3, Sec. 2.4) 

Receiver antenna height above the international ellipsoid 
(mean sea level) (earth radii) 

Relativistic Hamiltonian (Ref. P3, Sec. 3.2) 

Distance measured from the center of the earth to a 
field line (km) 

Parameter related to the longitudinal invariant (Ref. P3, 

Eq. (3.22)) 

Omnidirectional intensity of particles with energy greater 
than E (Ref. P3, Eq. (4.38)) 

Omnidirectional spectral intensity (Ref. P3, Eq. (4.32)) 
Modified Bessel function of order p 
Moment of inertia about the X' axis 
Moment of inertia about the Y' axis 


Pl-2 



Jj(x) 

j(c6,^>;0) 

J(O') 

J 

o 

i 

L 


M 


m 

m 

o 

N 

N(E) 


n. ,-n. 
l+l 1 


n 

P 


/s 

P 



P(f) 


r 

R 


Moment of inertia about the Z' axis 

First order Bessel function regular at the origins 

Unidirectional spectral intensity (Ref. P3, Eq. (4.32)) 

Directional flux density of particles as a function of 
angle between J and a plane perpendicular to B 

Omnidirectional flux density (particles cin sec" ) 

Length of nose cone and instrument section Explorer 4 

Angular momentum unit vector 

Mcllwain's geomagnetic shell parameter (earth radii) 

Earth dipole moment (Ref. P3, Eq. (1.9)) 

Relativistic mass of electron 
Rest mass of electron 

The number of digitizer counts per digitizing interval T 

Absolute differential spectral intensity, electrons/fission/MeV 
(Ref. P3, Eq. (1.5)) 

The number of digitizer counts measured to correspond to 
one frequency change (flip) on data record 

Index of refraction (Ref. P3, Eq. (4.10)) 

Relativistic particle momentum; also a transformation 
parameter in the conversion from B~L to R-\ (Ref. P3, 

Eq. (1.13)) 

Counter symmetry axis unit vector 

Component of particle momentum perpendicular to the 
magnetic flux density vector B 

Component of particle momentum parallel to the magnetic 
flux density vector B 

Relative spectral intensity of frequency f (Ref. PI, 

Eq. (Al.l)) 

Distance of a particle from the center of the earth (km) 

Vector from center of earth to tracking station in 
Vernal Equinox coordinates (km) 


Pl-5 



Vector from center of earth to satellite in Vernal 
l.quinox coordinates (Km) 

. 7 

Radius of the "arth, 0.63712 x 10 meters 

Equatorial distance from dipole center to B line (km) 

Extrapolated range of particles penetrating an absorber 
(Ref. PI, P3, P4; gm cm~2) 

Range for 50 percen^ transmission of radiation through 
an absorber (gm cm' ) 

Scaling factor - twice the number of accumulated particle 
counts corresponding to one frequency change on the data 
record (particle counts per flip-flop) 

Normalized particle momentum (Ref. P3, Sec. 4.23) 

Detector A response function (Ref. P4) 

Matrix element representing the kernel of the integral 
equation relating the true count-rate C(0) and the 
directional flux density -1(0') 

A parameter in the transformation from B-L to R-X 
(Ref. P3, Eq. (1.14)) 

Time interval selected for digitazation (usually 24 sec) 
Longitudinal drift period (Ref. P3, Eq. (3.2S)) 

Bounce period (Ref. P3, Sec. 3.26) (sec) 

Cyclotron period (Ref. P3, Eq. (3.9)) (sec) 

Bounce pe - *-’ id of a particle with mirror-point at 

latitude defined by Eq. (2.8) (Ref. P2) (sec) 

pro 

Time at which particle count-rate for interval between 
i^-jquency cnange (flip) is assigned (center of interval) 
(UT min) 

Initial time of time interval T (UT min) 

Vector from tracking station to satellite 
Alfven velocity (Ref. P3, Eq. (4.6)) (km sec - *) 

Speed of particle 

Total energy (Ref. P3, Sec. 3.2) 

Axis of maximum moment of inertia Explorer 4 



Axes of Vernal Equinox inertial coordinate system; X 
in direction of Ve_rnal Equinox of 1951; Y in earth's 
equatorial plane; Z to form right hand orthogonal system 

Axes of satellite fixed rotating coordinate system; X 
space fixed in direction of angular momentum vector L; 

Y y in the direction of satellite symmetry vector A 
(rotating in a_j>lane perpendicular to L in propeller¬ 
like motion); Z f to form a right-hand orthogonal system 

Axes of satellite inertial coordinate system; X in 
direction of angular momentum vector L; X -Y plane 
oriented to include B the magnetic flux lenlity unit vector 
at time zero when an r-f null is observed; Z to form a 
right-hand orthogonal system 

First order Bessel function irregular at the origin 

Atomic number of scattering atoms 

Axis of minimum moment of inertia Expioi-er 4 

Pitch angle. The angle between the particle momentum 
vector and the magnetic flux density vector (deg) 

Pitch angle of the particle that defines the loss cone 
(deg) 

Ratio of particle speed to speed of light in free space 

Relativistic mass ratio (Ref. P3), 

Spectral exponent (Ref. P3, Eq. (4.19)) 

Phase angle of tumble motion of satellite 

Geodetic latitude (radian) 

Geocentric latitude (radian) 

Geocentric latitude of satellite (radian) 


Phase angle of roll (spin) motion of satellite 


Detection efficiency of counter 


r d, 


(Aqh*\ 


Hour angle of the Vernal Equinox (radian) 

Angle between the Greenwich meridian and Vernal Equinox at t 
East longitude of the tracking station (radian) 



t Fuler an^le specifying the opening angle of precession 

about angular momentum vector L (Kef. PI, P4, Fig. 2,b) 


e 


e* 


\ 

H- 

r 


T. T 

L > p 

T 

l T> 


H 

K 

$ 

* 

$ 

*/> 


i/y(e') 


COp 

(*>L 



Angle between the directional counter symmetry axis and 
the plane perpendicular to the magnetic flux density B; 

0 = n/2 -a (deg); A = n/2 -0 

Angle between directional flux density vector J and the 
plane perpendicular to the magnetic flux density vector 
B (deg) 

East longitude of satellite fr^m Vernal Equinox (radian) 
Geographic latitude (radian) 

Conjugate geographic latitude (radian) 

Geomagnetic latitude (radian) 

Magnetic dipole moment 
Time parameter 

G-M counter deionization (deadtime) (sec) 

Scintillation counter deadtime 
Whistler frequency (Ref. P3, Eq. (4.13)) 

Cyclotron frequency (Ref. P3, Eq. (4.2)) 

Urn wave frequency (Ref. P3, Eq. (4.2)) 

Flux invariant (Ref. P3, Eq. (3.28)) 

Euler angle measured in X-Y plane 

Pi'ecession rate about angular momentum vector L (rad sec”*) 

Euler angle in the plane of rotation about Z'-axis from 
the line of nodes to the rotating satellite fixed X'-axis 
(deg) . 

Element of an orthonormal set 

Residual angular velocity about Z'-axis after transition 
to propeller-like configuration (rad sec” ) 

Plasma frequency (Ref. P3, Eq. (3.28)) 

Doppler shifted whistler frequency (Ref. P3, Eq. (4.7)) 

Angular velocity about Z'-axis (rad sec"*) 

Initial angular velocity about the Z'-axis at orbit in¬ 
sertion (rad sec ) 


Pi-6 
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Pig. P2-1. Block diagram. Explorer 4 Instrumentation and telemetry 
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A. Summary of Pn.D. Dissertation, Donald J. Manson, "Van Allen 

Radiation Belt and ARGUS Directional Flux Density Distributions, 
Explorer 4, Satellite Data 1 '; Saint Louis University, Physics 
Department (February 1967). 

Al. Introduction 

A recent determination (Ref. Al) of the correct roll (or spin) 
rate for satellite Explorer 4 has allowed, for the first time, system¬ 
atic and reliable determinations of directional flux density distribu¬ 
tions of geomagnetically trapped charged particles Injected by the 
ARGUS Event 1 of 27 August 1958. The results of the present analysis 
(1) indicate a pronounced and unique omnidirectional penetrating radia¬ 
tion effect associated only with ARGUS Event 1, (2) include a comparison 

-2 -1 -1 

of directional flux density distributions (particles cm sec ster ) 
of the naturally (inner V-n Allen Belt) and artificially (ARGUS Event 1) 
injected radiation, and (5) permit determination of a decay curve for a 
limited time period of the measured maximum (perpendicular to B) direc¬ 
tional particle flux density of the artificially injected (ARGUS Event 1) 
radiation. 

A2. Analysis and Experimental Details 

The angle 6 between the directional (scintillation Detector A, 
Channel 2,5) counter arid a plane perpendicular to the magnetic flux 
density vector B is determined as a function of time by calculating 
the rotational (roll and tumble motion) and forward motions of the 
satellite i.» the Vernal Equinox coordinate system. Critical phase 
angles for the roll and tumble motions at an initial reference time 



are determined by analysis of the r-f electric field strength patterns 
from the telemetry and by "loop closure-search" (based on the assumption 
that at equal angles of the counter axis relative to B the counting rates 
are equal; see Ref. A2 for details). Each short time-period section of 
count-rate data between signal fadeouts (corresponding to one-half a 
tumble cycle) is analysed separately. Satellite telemetry yields pri¬ 
marily count-rates versus time (Ref. A2) which, combined with the cal¬ 
culation of 0 as a function of time, yields observed count-rate C q ^ s as 
a function of 0. Correcting for deadtime, counter efficiency and effec¬ 
tive area (Ref. A2), and solving the integral equation (Ref. A3) relating 
count-rate to the directional flux and the counter response function, 

yields by an unfolding procedure the directional flux density J(0‘) 

-2 -1 -1 

(particles cm sec ster ). With the aid of the magnetic mirror con¬ 
dition J „(0'), (0', angle between vector J and a plane perpendicular 

ulclX 

to B), the directional particle flux density perpendicular to B is then 
calculated for larger values of B of the L shell on which the J(0') vs 
0' distribution data was acquired by Explorer 4. 

Figure A1 demonstrates the validity of the method used to calcu¬ 
late 0 as a function of time by indicating the consistency of all the 
data acquired by the directional counter for 4 minutes of trajectory 
in the Natural region of radiation (L = 1.35-1.44; approximately 1200 
miles of satellite trajectory). It is evident that much of the observed 
scatter is due to decrease in the maximum count-rate as the satellite 
moves to increased values of B and L, and that the shape of the count- 
rate vs 0 distribution is sensibly constant with small variations of L 
for the Natural radiation in this region. The relatively constant 
counting rate for large values of 0 indicates a penetrating component 
of the radiation. 
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A3. Results 


All the data acquired from the unshielded G-M counter (Detector C, 
Channel 5) and the directional scintillator (Detector A, Channel 2) for 
Passes 4l4 and 427 are summarized in Figs. A2 (Upper and Lower). Tri- 
angle-points for each of the two Figures are count-rate maxima (counter 
axis perpendicular to ll) determined from directional distributions such 
as shown by Fig. Al, dot-points are in turn the minima (counter axis 
parallel to 1'); the dashes are the average count-rate from the unshielded 
G-M counter (detects electrons of energy greater than 3 MeV; protons, 
greater than 31 MeV) and the X-points are the average count-rates from 
the shielded (Detector D, Channel 1) omnidirectional G-M counter (detects 
electrons of energy greater than 6 MeV; protons, greater than 40 MeV). 
Decrease in maximum count-rate with increasing B and L (Fig. A2) is 
evident for the Natural radiation (to the "left" of the arrows along 
the time axes); clearly visible is ARGUS Event 1 at L = 1.62 - L = 1.76. 
For Pass 4l4 a core of penetrating radiation is indicated (between arrows 
along time axes) at L = 1.71 to I. 76 , 1.8 hrs following the nuclear event; 
this directional penetrating component is .ot observed 24 hrs later {Pass 
427). Comparison of Passes 414 and 427 indicates essential constancy 
of the Natural radiation but rapid decay of the ARGUS shell within a 
24-hour period, and also indicates a change in energy distribution 
suggested by a difference in ratios (of the penetrating component of 
the directional scintillator count-rates to the average count-rates of 
the omnidirectional unshielded G-M counter) for the same 24-hour period. 

Figure A3 indicates a comparison of the directional flux density 
distribution, J(0 1 ) vs 4 \ for the Natural (left-hand plot), ~nd the 
ARGUS Event 1 radiation (right-hand plot) of 28 August 1958, 23.48 hrs 


following the event. 



A4. Short-lived Penetrating Radiation Phenomena Possibly Related to 
Geomagnet leally Trapped Debris 

This section contains a revised (relative to our Final Report, 
DASA Report l6l3, 15 January 1965 ) interpretation of the penetrating 
radiations observed during Pass 4l4 (over Huntsville) of Explorer 4. 
Hence Figs. A4, 5. 6 from Report l6l3 are reproduced in the present 
section. 

Figures A4, A5, and A6 are the result of detailed directional 
analysis of the data acquired from the region Indicated by the arrows 
on time axis (Fig. A2, Upper) where the penetrating radiation mentioned 
in A3 above was observed. The method of treating each limited portion 
of telemetry data between signal fadeouts allows plotting of the count- 
rate, vs time and 6 , during each fractional part of a rotation of the 
scintillator counter axis with respect to B. As indicated by Fig. A1 
this method yields a consistent disk-like shape for the Natural radia¬ 
tion with all data clustering around one curve. 

Figures A4, A5, and A6 indicate considerable departure of the 
ARGUS Event 1 radiation, 1.8 hrs after the event, from the disk-like 
shape of the Natural radiation and so each partial rotation of the 
counter axis with respect to"lT has been plotted separately except for 
Fig. A4 (b,c,d) where three such rotations (from the early penetration 
of the ARGUS shell have been superimposed to yield a count-rate vs 9 
distribution much like that for the Natural radiation). 

In Fig. A4(f) the time sequence shown by the arrows indicates 
Inal the ave -age count-rate is increasing rapidly during one partial 
rotation (approx. 3.6 secs) so that at equal angles of the counter 
axis with respect to B (e.g., 9 - oO deg) two different count-rates 



are observed, the later count-rate having increased as the satellite 
penetrated farther into the ARGUS shell. It should be noted that 
the omnidirectional penetrating component (at large 0 values) has 
increased rapidly (from 2000counts sec 1 in Fig. A4(b,c,d) to approxi¬ 
mately 5000 counts sec" 1 in Fig. A4(f).) The large count-rate (9000 
sec" 1 , 0 = 92 deg, Fig. A4(f)), at larger angles of 0 indicates a 
departure from the disk-like shape of the radiation. Figure A4(g) 
indicates that the scintillation counter is saturating (count-rate, 
greater than 10,000 sec" 1 ) when perpendicular to]? (0 = 0) while the 
count-rate at 0 = 50 deg has dropped (from 9000 to 7000 sec 1 ). 

Fir-re A4(h) indicates saturation at small to moderate 0 angles 
( > 10,000 sec" 1 , 0 40 deg) and a count-rate of almost 7000 sec 1 

at 0 = 60 degrees - again a considerable departure from disk-like 
shape of the count-rate vs 0 distribution. 

Figure A5(j) (acquired 2.5 secs later than Fig. A4(h)) indicates 
a more normal disk-like shape but for a limited range around 0 = 0 . 
Figures A5(j,k,l and ra) if superimposed would yield a single consistent 
distribution much like Fig. A5(m) where the C obg vs 0 distribution is 
disk-like (0 £ 30 deg) but with a large omnidirectional penetrating 
background component ( ^ 5000 sec 1 ). The final 5 points of the time 
sequence (arrows) shown in Fig. A5(m) indicate a sudden saturation of 
the counter (0 £ 30 deg). 

This saturation is continued in Fig. A 6 (n) at 0 = 30-35 degrees 
but is followed in Pigs. a6(o) by an abrupt change to a more nearly 
disk-like distribution; followed immediately in turn by saturation for 
all angles 0 ^ 40 deg. Figure a 6 (p) (0.2 sec later) indicates a 

return to the normal disk-like shape followed by a sharp drop in the 



count-rate at 6 - 0 as the satellite leaves U.a ARGUS shell. Figure 
At(q) is for the last measurable count-rates from the ARGUS Event 1 
wing as the satellite moves into the slot between the two Van Allen 
belts away from the ARGUS shell. 

The analysis presented relative to Figs. A4, 5. 6 suggests that 
the data from Pass 414, 1.8 hrs after Event 1 contains evidence for 
short term penetrating components (directional and omnidirectional) 
of the radiation which are neither spatially homogenous nor Isotropic. 
The fact that the scintillation counter goes to saturation suddenly 
during one portion of the rotation of the counter axis with respect 
to vector B as in Fig. A5(m) and A6(o) while the remaining portion 
of the rotation measures a fairly normal distribution suggests a vary¬ 
ing directional component of the flux at these points (possibly due to 
gamma radiation from trapped debris near the ARGUS shell). The normal 
disk-like distributions (Fig. A3) measured in the ARGUS shell 1 day 
later for Pass 427 apparently are present in Pass 4l4 but are obscured 
and distorted by the large penetrating omnidirectional component and 
bursts of directional radiation which saturate the counter. These 
effects are greatly diminished in Pass 415, (3.67 hrs following Event 1) 
and are not found in Pass 427 (25.48 hrs following Event 1). 

No such effects were found from ARGUS Event 2 where data from 
Pass 454 (I .87 hrs following Event 2, 30 August 1938; Refs. A2, A4) 
indicate the usual disk-like distributions. 

Figure A7 indicates the time decay of maximum directional flux 
density f° r two d&y s following ARGUS Event 1. It was not 

possible to follow the decay for a longer period than »,wo days because 
by the third day tne directional trapped radiation intensity has dimin¬ 
ished to near background value. 

Ar, 



The x-points lie well on a straight line (exponential decay law), 
the solid-line curve is a plot of the t law (observed for the 

omnidirectional radiation, see Part C of the present report) for com¬ 
parison. Because of the snort time Interval represented, about 2 
days, a decision between the two decay laws is not considered valid. 

Figure A8 is a 3-dimensional plot of J (0'), the directional 
° max 

flux density perpendicular to vs B and L; both Natural and ARGUS 
Events 1 and 2 radiations are shown. Figure A8 summarizes most of 
the data analysed for the scintillation counter, Detector A. 

A5. Conclusions 

Detailed analysis of ARGUS directional flux density distributions 
from Explorer 4 data (Passes 414 and 427, over Huntsville, Alabama, on 
27 and 28 August 1938 respectively) for 1.8 and 25.48 hrs respectively 
after Event 1 and comparison with Natural flux density distributions 
indicate that the disk-like shape as measured by the scintillation 
counter in the ARGUS shell is not substantially different from the 
Natural radiation (by Fig. A3, the half-width at half-maximum fhwhm) 
is 7 deg broader on the average for the ARGUS J(0') compared to the 
Natural radiation; maximum values of directional flux densities are 
12 x 10 7 for the Natural vs 6.7*10 7 cm ^sec ^ster"* for ARGUS). The 
27 August 1958 data from Pass 4l4 (Ref. A2) and Pass 415 (Ref. A3) 
indicate that these measured distributions are temporarily obscured 
and distorted by a large penetrating omnidirectional component and 
a varying directional component. The short-time limited decay data 
suggest an exponential decay of maximum directional flux density at 
L = 1.7 but a t law observed for the omnidirectional data from 



G-M counter (Detector C) Is not excluded. It rj«ay be suggested 


that particle and energy differences between Natural and ARGUS radia¬ 
tions can be expected to yield somewhat different decay laws. 
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Pig. Al. Pass 427 (over Huntsville, Alabama, 28 August 1958) 

observed directional count-rate C vs 0 for 4 min 

obs 

in the region of Natural radiation. 

o, C obs vs 0; 227-5 - 229-5 min UT; L, 1.35 - 1-42 
earth radii; B, 0.149 - 0.164 gauss; 4 , C Qbs vs 0; 
229.5 - 231.5 min UT; L, 1.42 - 1.49 earth radii; 

B, 0.164 - .180 gauss. 

Data were plotted by combining about one-half of 
the 100 individual sections of C Qbs vs 0 described in 
Ref. A2, Sec. 3 . 3 . All data acquired within the time 
period 227.5 - 231.5 UT are plotted to show the amount 
of scatter in the data and to demonstrate the consist¬ 
ency of the analytic interpretation using the experi¬ 
mentally determined parameters: roll and tumble phase 
angles 8 and V , and the roll and tumble rates. 
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Pig. A2. (Upper, Lower) 

Natural and ARGUS Event 1 count-rates vs time and 
B-L coordinates for Passes 4l4 and 427 over Huntsville, 
Alabama. Comparison of Average Omnidirectional Count- 
Rate with Maximum (perpendicular to’B) and Minimum 
(parallel to"£) Directional Counting Rates. 

Triangle-points, maxima; and dot-points, minima, 
scintillation counter (Derector A) observed count-rates; 

Dash-points, 15-sec averages of observed omnidirectional 
(unshielded G-M counter. Detector C) count-rates; 

X-points, 6-sec averages of observed omnidirectional 
count-rates (shielded G-M counter, Detector D). 

Arrows on the time axis indicate time region in which 
distributions of Pigs, A4, 5. and 6 (Pass 4l4) were acquired 
(ARGUS Event 1). Calibration information for the three 
counters is summarized in the Preface (PI) of the present 
report. 
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Pig. A3. Directional flux density distributions J(0') vs 0 
for Pass 427, 28 August 1958; Natural (left-hand 
plot) and ARGUS Event 1 (right-hand plot) radiations. 
Variations of magnitudes of «J max (0') = J(0) and shape 
of the distributions are shown. J(0') for ARGUS 
Event 1 radiation were corrected for varying penetrating 
radiation component as explained in Ref. A2. 
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Pigs. A4, A5, and A6 

Detailed sequence of scintillation counter count-rates 
Cq^s vs ® (angle between scintillation counter axis 
and the plane perpendicular to *B) for Pass 414 ARGUS 
Event 1 penetration, 27 August 1958* Where used 
arrows indicate time sequence in acquisition of points. 
Distributions indicate a large omnidirectional pene¬ 
trating component and a varying directional component 
of the radiation (discussed in detail in the text). 
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ARGUS Event 1 






ARGUS Event 1 









ARGUS Event 1 





Pig. A7. Decay of maximum directional particle flux densities 

J (0*} vs time for ARGUS Event 1. Plotted points 
max 

are for flux density distributions J(0*) at 0'= 0. 

The maximum ARGUS Event 1 radiation is observed to 
decay one order of magnitude in 45 hrs. The lower 
curve is a plot of the t law (observed for the 

omnidirectio-'1 radiation, Part C of the present 
report, for the same ARGUS Event 1 radiation). 
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Pig. A8. Directional maximum particle flux density J max ( 01 ) 
vs B and L for both Natural and ARGUS events 1 and 
2 radiations (Detector A, Channel 2 ) 


Detailed information is tabulated as follows. 


Natural Radiation 


L Slope Pass No. Date 

(earth radii) 


oL 

p 


1.360 

r . ... 

- 7.92 


28 Aug. 

1938 


1.380 

-4.82 


27 Aug. 

1938 



1.407 

-7.15 

427 

28 Aug. 

1958 

o 


1.409 

- 2.76 

414 

27 Aug. 

1958 

£ 


1.434 

-7.13 

454 

30 Aug. 

1958 

& 


1.458 

-7.19 

454 

30 Aug. 

1958 

7 


1.580 

-7.92 

427 

28 Aug. 

1958 


ARGUS EVENT 1 


e 


1.690 

-3.78 

414 

27 Aug. 1958 

Z 


1.694 

-6.42 

427 

28 Aug. 1958 

K 


1.724 

-4.10 

427 

28 Aug. 1958 

F 


1.728 

-4.11 

440 

29 Aug. 1958 

■? 


1.730 

-3.58 

415 

27 Aug. 1958 

5 


1.743 

-5.05 

427 

28 Aug. 1958 

0 


1.748 

-2.00 

427 

28 Aug. 1958 

X 


1.777 

-4.48 

427 

28 Aug. 1958 


ARGUS EVENT 2 


IT 


2.124 

-7.55 

454 

30 Aug. 1958 

P 


2.124 

- 3.60 

454 

30 Aug. 1958 

cr 1 


2.144 

-15.40 

454 

30 Aug. 1958 

r 


2.172 

-10.40 

454 

30 Aug. 1958 
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B. Summary of M.3. Thesis of Joseph F. Fennell, "Van Allen Belt and 


ARGUS Directional Flux Density Distributions, Explorer 4 Satellite 
Data, ARGUS Event 2", Saint Louis University, Physics Department 
(July 1966 ). 

Bl. Introduction 

This Thesis, Reference Bl, contains essentially an analysis of the 

-2 -1 -1 

directional distribution of flux densities (particles cm sec ster ) 
vs the angle between the counter symmetry axis and the plane perpendi¬ 
cular to the magnetic flux density vector, of geomagnetically trapped 
naturally occurring and artificially injected charged particles asso¬ 
ciated with ARGUS Event 2 for Pass 4^4 over Huntsville, Alabama (for a 
time interval of about 10 min. beginning 1.8 hrs following the Event). 

The purpose of this limited analysis was to compare the situation with 
that of Pass 4l4 (Manson Dissertation, Ref. B2) with the Passes occur¬ 
ring approximately 3 days apart. The count-rate data were obtained 
from a scintillation photomultiplier detector (Detector A, Channel 2 
Ref. B2); window aperture, 15 degrees half-width at half-maximum; 

A g ff£, 0026 cm ; deadtime 100 microsecs aboard Satellite Explorer 4. 

Tne thesis represents the first detailed examination of the direc¬ 
tional data for ARGUS Event 2 by the Saint Louis University group. There 
is presented a short discussion of a method for determining the proper 
roll phase angle used in this analysis and the results from determinations 
of satellite roll frequencies for times at which they had not been pre¬ 
viously determined with sufficient precisio-. to yield reliable determi¬ 
nations of observed count-rate vs an 6 ^- e ® (between the counter 

symmetry axis and a plane perpendicular to the magnetic flux density 
vector B). 

Bl 



B2. Determination of the Phase Angle for the Roll Motion of Satellite 


To determine the directional flux density J(@') vs angle 0' between 

the counter axis and the incident flux direction, C . (0) as a function 

obs 

of 0 must be known. 9 at time t (relative to an arbitrary reference time) 
is given by 

6 - A rc s i n[ b, S inOpi< *8) - b x Cos Cpt + Sj S incfit +Y ) 


+ b 3 Cos +S) Cos(4>t +Y)] , 


(Bl) 


where all parameters are known explicitly except 3 and 'Y , the 
roll and tumble motion phase angles. W was determined as 

7 - Arccos (A ‘ Ys), w) 

\ 

where A is the satellite symmetry axis vector and Y s is a unit vector 

in Vernal Equinox Coordinates (Ref. B2). 

3 is determined (Ref. Bl) in Eq. (Bl) by assuming counting 

rates are equal approximately at equal angles 0 in the Natural radiation 

region (over short periods of time, t < 0.5 min). Taking two times 

t , t. with equal counting rates, it can be shown (Ref. Bl) from Eq. (Bl) 
* J 

that 

S|j - Arcta.n ("fjj/j/j) (B3 > 


where f, 4 = f, + f., g,, = g, -*• 6.> 

1 J . 

f = b (t ) SinCl^t^ + P(t 1 ,y ) Cos{ Ip t i ), 

= b 1 u i ) cosn/ztj - ) sm{ \p t i ) 

y ) = b (t 1 ) Cos( <p> t i + y )-b 2 (t jl ) Sin( <p> t i+ y ). 


and 
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An average <j is then determined by considering 

§a.v<] = A ret a_n (" L (t/j /9ij )) 

N 


(B4) 


where the sum is over N time pairs t^, tj. The consistency of the data 
can be seen by considering the variations in the 3 ^ from 3 avg. 

If the spread in the values of 3 .. is small the random fluctuations 

3 

in the data are small. ^ivg was determined (with respect to an arbi¬ 
trary reference time) as 94° for Pass 454 and this value reproduced the 
count-rate variations in a simulator program (see Ref. B2). 

The roll rate for Explorer 4 was determined for the entire tele¬ 
metry times of Passes 4l4 through 454 by fitting a least square poly¬ 
nomial to the known roll rate values at discrete reference times of 
Passes 414, 415, 427, and 454. A second method for determining the roll 

rate, based on the fact that Sin 8 is a minimum when C . (8) is maximum 

obs 

is presented by Fennell (Ref. Bl). (The second method was applied to 
Pass 440 and the resulting roll rate was used by J. M. Paikeday (Ref. B3) 
to obtain a directional flux density distribution for Pass 440 ARGUS 
Event 1 radiation.) 


B3. Results 

With 8 determined as a function of time, plots C Qbs (8) vs 6 are 
made and then corrected for the counter deadtime and (effective area 
x counter efficiency) to yield the true count-rate C(8) 


C (9) = 


Co bS (&) 


(B5) 


eAeff (1 r vCobs(d)) 



The corrected or true count-rates C( 0 ) are used to compute the corres¬ 
ponding directional flux densities J (9 1 ) as described by J. M. Paikeday 
(Ref. B3). Representative examples of the observed count-rates C Q ^ s (0) 
vs 0, true count-rates C(0) vs 0, and directional flux densities J(0 1 ) 
vs 0', for Pass 454 Natural radiation are shown by Pigs. B1 and B2. The 
C o b S (0) vs 0 curves showed an "open loop" phenomena (see Pig. B3(a)) 
which is due in part to the rapid rate of change of both the omnidirec¬ 
tional and directional components of the radiation as the satellite moves 
through the ARGUS Event 2 shell. The curve was reduced to a single valued 
curve by subtracting out this component (compare Pig. B3(a) and 
of Pig. B3(c)). The resulting C Q ^ s (0) vs 0, C(0) and J(© 1 ) vs 0, 0 1 
curves for the ARGUS radiation are shown in Figs. 33,4. Figures Bl,2,3,4 
and Tables Bl, B2 in the Pigs. B2, B4 captions represent the results of 
this analysis. 

B4. Conclusions 

The maximum directional flux density and observed count-rate vary 
considerably in the six minutes of the Natural radiation observed for 
Pass 454 while the half-width at half-maximum (hwhra) remains relatively 
constant at 17 + 5 degrees (see Table Bl, Pig. B2 caption). The varia¬ 
tion in maximum count-rate and flux density is probably related to the 
large spatial and magnetic flux density changes involved in the satellite 
motion. 

The non-zero J(0*) for 0' 50 bo 60 degrees (Pig. B2), in the 

Natural radiation region is due probably to pene rating (independent of 
0 1 ) radiation capable of transmission by the counter shielding (E* 3 Mev); 

it is the background count-rate for the detector in the region of Natural 
radiation. 
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The flux densities and hwhm vary considerably for the ARGUS Event 2 
Pass 454 radiation. By Pig. B4 the ARGUS radiation appears generally to 
have a smaller hwhm ( 16 ° vs l8°) and a larger maximum flux density 
(6.4 x 10 7 vs 4.9 x 10 7 particles cm ^sec -1 ster"^) than the Pass 4^4 
Natural radiation. It is noted that there are large differences in B 
and L coordinates between the Natural and ARGUS radiations. The results 
for Pass 454 are included in Table B2 for comparison with previous analy¬ 
sis (Refs. B2, B5, B4, 35). 

The maximum flux intensity and hwhm of the Pass 454 ARGUS Event 2 
are considerably lower (6.4 x 10' vs 2.3 x 10° particles ctn sec^ster -1 ) 
and narrower (16° vs 2lf°) respectively than those of Pasa 4l4 for ARGUS 
Event 1 with both Passes occurring about 1.8 hours after their respective 
events. Assuming both nuclear devices were the same, the differences in 
Pass 414 and Pass 454 ARGUS results may be due to spatial or geomagnetic 
affects. The data occurred at the two L's, I .72 (Event 1) and 2.1 
(Event 2). 

Some tentative conclusions can be offered for Pass 454 Event 2, as 
follows: 

(1) ARGUS radiation is more intense than the Natural as is evidenced from 
the Tables B1 and B2 (see Pigs. Bl, B4 captions). 

(2) ARGUS radiation appears to be slightly more disk-like than Natural 
as is evidenced from the hwhm (Tables Bl, B2, B3). 

.3) ARGUS radiation showed no sign of any rapidly varying omnidirectional 
penetrating radiation effect as was observed by Manson (Ref. B2) in 
Pass 4)4 
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Since Pass 4^4 is the only ARGUS Event 2 directional data which 
has been analyzed by the methods presented In Refs. B1 and B2 to date, 
a more intensive study of Event 2 data and a more complete comparison 
with Event 1 directional data must await the analysis of more Event 2 
directional data. 
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Pig. Bl. Representative ° ol3S (®) vs ® curves for Pass 4^4 Natural 

radiation. Actual data points are shown to indicate the 

scatter in C . (6) vs 9. 
obs 
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Fig. B2. True count-rate C(0) (counts cm sec corrected for deadtime, 

effective area and counter efficiency) and directional flux 

-2 -1 -1 

J(0') (particles cm sec ster ) vs 0 and 0' for the Fig. B1 
C obs (e) curves. The bars on the J(0') vs 0* curves indicate 
the scatter in calculated values of J(0*) due in turn to the 
scatter of points in the corresponding C Qbs (0) vs 9 curves. 


Table Bl. Characteristics of Natural radiation (from Ref. Bl) 


Time (min) 

*hwhm (deg) 

J(6') 

max 

B (gauss) 

L( earth radii 

300.00 

22.0 

4.9x10^ 

.1734 

1.376 

301.221 

17.5 

4.9x10^ 

,1782 

1.414 

501.913 

17.5 

4.7x10; 

.1806 

1.434 

302.370 

17.5 

3.5xlo; 

.1834 

1.458 

304.570 

16.0 

1.4x105 

.1912 

1.583 

305.725 

16.0 

8 .2x10^ 

.2172 

1.640 

306.990 

12.0 

7-7x10 

.2316 

1.776 


* half-width at half-maximum 
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Pig. B3. 

( a ) C Qbs (e) vs 9 ARGUS Event 2 'loop' before the linearly 
increasing background has been subtracted. Compare 
Pig. B3(c) which is the resulting C Qbs (0) vs 9 curve 
after subtraction of background. The background sub¬ 
tracted is represented by (constant x t,), where t^ is 

the time (following the time of the first point of a 

t/h 

given C Qbs (0) vs 6 curve) corresponding to the i° n 
point of the given curve. 

(b)(c)(d) c 0 b S (®) vs 6 f° r Pass ^5^ ARGUS Event 2 radiation 

after loop closure is accomplished. 
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Pig. B4. C(0) and J(0') vs 0 and 0’ for the C Qbs (e) vs 0 curves 
(Pig. B3), Pass 4^4 ARGUS Event 2. The bars on curves 
indicate the scatter in calculated values of J(0') as 
above (Fig. B2). 

Table B2. Characteristics of the Pass 4^4 ARGUS Event 2 radiation 
(taken from Ref. Bl). 

Time (min) hwhm (deg) J(0') v B (gauss) L(earth radii) 

fTlclX 


310.178 

12.0 

3.2xl0 7 

.2720 

2.120 

310.223 

24.0 

5.65X10 7 

.2726 

2.124 

310.279 

13.0 

5.95X10 7 

.2740 

2.128 

310.336 

14.0 

6.45x10 7 

.2742 

2.144 

310.496 

19.0 

2.82X10 6 

.2764 

2.172 
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C. Summary of Ph.D. Dissertation, John A. George, “Omnidirectional 
Fluxes; Explorer 4 Satellite Data, ARGUS Events 1 and 2”, Saint 
Louis University, Physics Department (February 1967). 

Cl. Introduction 

In August and September 1958, three low yield nuclear detonations 
(ARGUS Events 1, 2, 3) were produced at high altitude over the south 
Atlantic Ocean. Electrons emitted were subsequently trapped in the 
geomagnetic field. The characteristics of this radiation were measured 
by instruments on board Jason sounding rockets and in the Satellite 
Explorer 4. A number of papers have appeared in which the results of 
what is known as the ARGUS Experiment were described (Refs. Cl, C2, C3). 

The present Dissertation is concerned with further analysis of 
data obtained by Explorer 4 for ARGUS Events 1 and 2 for the G-M 
Detectors C (Channel 3) and D (Channel 1); therefore data for omni¬ 
directional radiation is considered. 

C2. Experimental Details 

The characteristics of satellite Explorer 4 (1958 Epsilon) have 
been described in detail by Van Allen (Ref. C2). For later calibration 
information, as used in the present report, see Preface P2 of this 
report. Hie Anton 302 G-M counters only are involved in the present 
analysis. A summary of the characteristics of these two detectors is 
given in Table Cl. These detectors are omnidirectional and do not 
differentiate the kinds of particles being detected. However for the 
ARGUS radiation, it will be assumed that the predominant contribution 
to the counting rates is b> electrons resulting from the fission 
processes. Anton 302 G-M counters similar to those on Explorer 4 


Cl 



have since been flown on a number of satellites, and it has been con¬ 
cluded (Refs. C4, C5, C6) that the counting rates must be attributed 
to penetrating electrons of E e > 1.5 MeV. Motz and Carter (Ref. C7) 

studied the response of Anton 302 counters to fission-spectrum 

2 

electrons, They found that up to 2 gm/cm of lead-absorber thickness 

the curves were as expected for 6-transmission; however, for energies 

beyond this range the transmission was enhanced somewhat by bremm- 

strahlung. Detector C (the "unshielded" G-M counter, Channel 3) is 

2 

shielded by 1.2 gm/cm of stainless steel, which is equivalent to a shleld- 

2 2 
ing of 1.46 gm/cm of aluminum. The extrapolated range Rg (gm/cm ) 

versus energy E(MeV) for aluminum given by Marshall and Ward (Ref. C8) 

is 


R,. = 0.528E -0.094 (Cl) 

For Detector C, Rg by Eq. (1) corresponds to 2.95 MeV electrons 

(Table C2); for Detector D (the "shielded" G-M counter, Channel 1) 

2 

which had an additional shielding of 1.6 gm/cm of lead, the total 

2 

shielding is equivalent to 3.89 gm/cm of aluminum and Rg corresponds 
to 7.6 MeV electrons. For the unshielded Detector C at least the 
assumption of electrons being counteJ directly seems quite valid. How¬ 
ever, the count-rate of the shielded Detector D may be enhanced some¬ 
what by bremmstrahlung. A plot of the omnidirectional geometrical 
factors versus the energy necessary to penetrate the various thickness 
of material, using the information of Tables Cl an< C2, is given in 
Fig. Cl. 

An automatic digitizer was used to analyze the telemetry play¬ 
back data and transform it to punch cards. A computer program was 
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written to obtain count-rate versus time for each selected Pass of the 
Satellite fxplorcr 4 over Huntsville, Alabama. The observed count- 
rates C Q k s (sec ) were corrected for the deadtime of the detectors and 
associated electronics to obtain true count-rates C (sec j by using 
the approximation (instead of hq. (Dl), p. U2 below) 



where deadtime *£« is 62.5 microsecs (Ref. C2) and was confirmed in 
the present work. Figure C2 is typical of the data available for 
analysis. 

The B and L coordinates for the ARGUS region were obtained from 
GE TEMPO Data Center (Santa Barbara, California) through the offices 
of the Defense Atomic Support Agency. Figure C3 shows the regions of 
interest traversed by Explorer 4 in B-L space for typical Passes over 
Huntsville, Alabama. Table C3 lists the complete set of Passes 
analyzed and are arranged in three columns, each of which corresponds 
to nearly the same region of space relative to the earth. Tne orbi¬ 
tal frequency of Explorer 4 was about 13 orbits per day. 

In order to investigate the possible correlation of Explorer 4 


ARGUS count-rates with solar activity. 


3-hr K 

P 


indices are presented 


by Fig. C4 for the period under investigation 27 August to 4 September 


19S8. Subsequent to the large geomagnetic fluctuations observed at 


the onset, a period of decreasing activity to relative magnetic quiet 


levels followed. No sudden commencements were observed until 5-5 


September, a period which is classified as one of magnetically dis¬ 
turbed days. 





C3. Results 


ARGUS Event 1 

The decay of the fission produced electrons trapped in the earth's 
magnetic field is indicated in Fig. C5, a log-log plot of the true 
count-rate versus T , the elapsed time since Event 1, at a fixed L 
value. The plotted points are the true omnidirectional count-rates 
of the unshielded G-M counter (Detector C, Channel 3). The count-rates 
are differences between total and estimated background rates, and are 
for the 414, 427, ...518 sequence of Passes. The maximum intensities 
observed are well represented by the L = 1.72 shell. The peak location 
did not shift, at least within the experimental error, from this L 
shell during the time interval indicated. If it is assumed C (sec *) 
versus t(hrs) is given by 

C = C Q t“ n , (C3) 

where C q and n are undetermined constants, a least squares fit of the 
data for L = 1.72 yields n = 1.13, C q = 1.31 x 10 3 ; whence Eq. (C3) 
becomes 

C = 1.31 x 10 5 t' 1,13 , 0.19 < B <0.23. (C4) 

This t~*'* 3 law agrees quite well with the decay rate determined by 
Cladis and Walt (Ref. C3) from Jason rocket data for ARGUS Event 2 
(L = 2.1). 

Figure C6 is a plot of the count-rate ratio unshielded/shielded 
omnidirectional G-M counters. Data for shielded and unshielded G-M 
counter of rates greater than 10 sec * and 30 sec ^ respectively, only 
were used to eliminate data averaged over larger regions of space. 
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The data for Pass 414 may be in considerable error due to saturation 
of the counters. The vertical bars are estimates of the spread of data 
arising principally from the scatter in the count-rates of the shielded 
G-M counter. The count-rate ratios terminate at Pass 466 due to the 
difficulty of separating the ARGUS electrons from the natural back- 
ground radiation. The dashed line is the plot of C = C Q t * . Of 

interest is the sudden broadening of the ARGUS Event 1 shell subsequent 
Pass 453, Fig. C7. The shell width at half-maximum is given in terms 
of ^ L.Pass 453 was the first pass *0 record the effect of Event 2 
but all data of Fig. C7 are for Event 1. However, Event 2 occurred 
subsequent to the satellite traversing the Event 1 Pass 453 shell. 

Since Passes up to 492 occurred during a period of low geomagnetic 
activity (K^ < 3 q ) it is suggested that the broadening observed on 
Fig. C7 is due to ARGUS Event 2. 

Figure C14, a 3-dimensional plot of true count-rate vs B and L, 
summarizes the data analyzed for the unshielded G-M counter (Detector C, 
Channel 3) for ARGUS Event 1. 

ARGUS Event 2 

The data analyzed for ARGUS Event 2 is over a smaller time inter¬ 
val than Event 1. Figure C8 gives the true count-rate versus elapsed 
time since Event 2. Again the data is for count-rates above the 
estimated background. The peak intensities were centered veiy nearly at 
L = 2.11. Pass 453 was the first to indicate ARGUS Event 2, 0.03 hours 
after the event. Considering the location of the satellite (297° E, 

27° N) and the detonation (352° E, 50° S), fission electrons would need 
to have energies in excess of 19 MeV to reach the satellite position in 



their eastward drift in the 108 seconds (Ref. C25), Although there 
is a significant numoer of MeV-range 8-electrons available from the 
8-decay spectra of fission products, the 8-decay end point energies are 
well below 19 MeV. Hence the initial count-rate may be due to direct 
illumination by 'y -rays emitted by the fission debris. The bomb 
debris would have to expand to a higher altitude than the detonation 
point (about 480 km) in order to be in the line of sight of the satel¬ 
lite. Colgate (Ref. C9) discusses the mass motion of a high altitude 
nuclear explosion. He assumes that the detonation products produce a 
cavity in the magnetic field and, in general, the hole will expand 
until the work done against the magnetic field equals the kinetic energy 
of the detonation products. When the expansion of the detonation plasma 
products ceases, the subsequent motion will be a spreading along the 
lines of force away from the detonation hemisphere. This occurs in the 
order of seconds according to Colgate. With the debris distributed 
along the lines of force for L = 2.11, the satellite would be in a 
position to detect the emitted y -rays. 

The decay rates of ARGUS Event 2 are very nearly represented by 
a 1/t law if the data for Passes 466-492 are considered. However, 

Passes 505 and 518 indicate increases in count-rate in a manner diffi¬ 
cult to explain (Fig. C9). A sudden commencement of solar activity 
did occur shortly after Pass 505 and 3-5 September are listed as mag¬ 
netically disturbed days. 

It is noted that the shielded G-M counter (Detector D, Channel 1) 
had very low count-rates for the sequence of Passes analyzed for ARGUS 
Event 2. These were typically less than 20 sec*. Figure CIO shows 
the count-rate ratio of the unshielded/shielded G-M counter for the L 
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shell of peak intensity (L = 2.11) for Passes 466-492. Subsequent to 
Pass 492 the ARGUS Event 2 electrons detected by the shielded counter 
were indistinguishable from the natural radiation. 

C4. ARGUS Loss Mechanisms 

The effect of the atmosphere on the decay of the trapped electrons 
has been calculated using a Fokker-Planck Reproach to derive a diffusion 
equation (Ref. CIO). This technique was applied by :.*alt (Ref. Cll), also 
Walt and Newkirk (Ref. C12) to the observed time variation produced by 
the Starfish nuclear explosion. In regions of the magnetosphere for 
which L = 1.3, it was found that the dominant factor limiting the 
trapping time of electrons is scattering by components of the atmosphere. 
However, for L = 1.3 the decay of trapped electrons is dominated by 
other loss mechanisms of electromagnetic-hydromagnetic character. 

A breakdown of the adiabatic invariants which will eventually 
enhance the atmosphere losses has been proposed by a number of authors. 
The processes responsible for the violation of the third adiabatic 
invariant have been considered by Parker (Ref. C13), Davis and Chang 
(Ref. C14), and Nakada and Mead (Ref. C15). Large scale geomagnetic 
fluctuations associated with sudden commencements and the initial phase 
of geomagnetic storms can cause diffusion of electrons across L shells. 
However, the dominant effect would be in the outer radiation zone, 
and repeated magnetic storms would seem to be necessary for rapid dif¬ 
fusion. The exact analytic details are quite complicated and have not 
been attempted so far in the present work - however the concept is 
highly suggestive. Sudden commencements were observed on August 27 
at 0243, 0303, and 0541 hrs UT. Pass 414 traversed the ARGUS Event 2 
region at approximately 0456 UT, 1.8 hours after the event. No further 
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sudden commencements were observed until September 3 and 4 which co¬ 
incide with Passes 505 and 518 respectively. 

The effect of several other electromagnetic-type loss contribu¬ 
ting factors have been studied by workers in this field, as follows. 
Resonant particle-wave interaction in which the magnetic moment in¬ 
variance is violated has been studied by Wentzel (Ref. C16), Dragt 
(Ref. C17), and also Chang and Pearlstein (Ref. C18). A summary of 
Dragt‘s work (Ref. C-17) on the effect of hydromagnetic waves is given 
by Fig. Cll. Critical energies below which the magnetic moment is con¬ 
served in interactions with hm waves are shown for particles mirroring 
at latitudes greater than 30 degrees assuming the hm wave frequency is 
cut off sharply at v* mgx » Various cut-off frequencies are indicated 
as it is assumed that hm waves may vary from 1 to S sec * up to the 
local proton cyclotron frequency. It is clear that such frequencies 
would not be very effective in interaction with ARGUS electrons. 

Kentzel gives a qualitative discussion to the effect that hm waves could 
influence the mirror point distribution of the relativistic ARGUS 
electrons during a geomagnetic storm. This may be a contributing 
factor to the large loss observed within the first day which was active 
magnetically. 

Again, the precipitation of electrons through the action of 
whistlers has been investigated by Dungey (Ref. C19) and by Cornwall 
(Ref. C20). For resonant interaction to take place, the Doppler- 
shifted frequency of the whistler as seen by the electron must be equal 
to the relativistic cyclotron frequency of the moving electron. Figure 
C12 summarizes Cornwall's results showing the L dependence of the 
precipitating particles with increasing electron energy. The ARGUS 
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Events 1 and 2 regions a^e shown, indicating, at least qualitatively, 
that whistlers ir.ay not be a major mechanism driving electrons to 
larger B values. 

Also, Faulikas et al (Ref. C21) suggest longi'.udinal waves 
propagating along magnetic field lines would interact with electrons at 
their mirror points and would be an efficient loss mechanism. These 
electrostatic waves may be able to move electrons to larger B values 
in times of the order of days. 

Further detailed theoretical investigations are needed before 
quantitative results may be given to adequately explain the rapid 
decay of the ARGUS electron shells. 

C5. Conclusions 

The analysis of the data obtained by the omnidirectional G-M 

counters (Detectors D, C; Channels 1, 3) carried by Explorer 4 for the 

ARGUS Events 1 and 2 can be summarized as follows. 

The two Events are centered on the L = 1.72 and L = 2.11 shells 

respectively. The locations of the peak intensities did not vary, at 

least within experimental error, from these values during the time of 

-1 13 

observation. The ARGUS Event 1 region showed a decay law of t 

for E > 3 MeV. The count-rate ratio of the unshielded to the shielded 
e 

counter was found to decrease with time. This indicates spectral 
changes occur which is in apparent contradiction with the results of 
the published Jason rocket data. A general broadening of the ARGUS 
Event 1 region was observed subsequent to the detonation of ARGUS 
Event 2. The increase in the width at half-maximum was by a factor of 
1.5. ARGUS Event 2 seemed to be influenced to a greater degree 
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than ARGUS Event 1 by the magnetically disturbed days of September 3-S. 
Mcllwain (Ref. C22) theorized that magnetic storms produce a decrease 
in electron flux which occurs so quickly that it is difficult to follow 
and then is followed by an increase larger than normal. The observed 
increases on Passes 505 and 518 of Event 2 may be related to these 
magnetic storm effects. The larger value L shells of Event 2 compared 
to Event 1 may account for the variation of the storm influence. Al¬ 
though a general decrease in count-rate ratios were observed for both 
events, ARGUS Event 2 was not as sharp as the first event and may point 
to L dependence of the decaying fission spectrum. The data is too 
limited to determine the precise nature of this dependence. 

The decay of the ARGUS Event 1 region was found to follow a 

-1 13 

t ' law. It may be coincidental that this is nearly identical to 

the fission product decay law. The fission fragment 0-ray source 

-1 2 

strength decays approximately according to a t ' law (Ref. C23). 

-1 2 

Measurements indicate that although the t * decay represents a 

reasonable average there have been instances where exponents in the 

range -0.9 to -2.0 are required to represent the rate of decay. 

Figure 13, from Heller (Ref. C24), gives the number of beta rays per 
235 

fission from U per energy interval versus time for various energies. 
The ratio of high energy to low energy y3-rays decreases with time. 

This contradicts the results of Figs. C6 and CIO which indicated the 
unshielded/shielded count-rate ratio that is, C(E e >3 MeV)/ 

C(E g > 7.6 MeV). decreases with time. Further study is required on 
the lifetimes of the fission debris which is trapped in the earth's 
geomagnetic field. 

A brief study was conducted of possible mechanisms which would 
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explain the rapid decay of the ARGUS Events. Theories involving the 
breakdown of the adiabatic invariants have been proposed by a number 
of workers. The author (J. George) applied these theories to the 
ARGUS data. Computations indicate that hydromagnetic waves and whistler 
interactions resulting in the violation of the first invariant are not, 
at least qualitatively, major mechanisms in removing ARGUS electrons. 
Although resonant particle-wave interactions, and other processes such 
as diffusion, may operate, they do not dominate. Further theoretical 
work is necessary to determine the controlling mechanism for the ARGUS 
decay. 

Adiabatic processes such as betatron acceleration and decelera¬ 
tion cannot be excluded in any detailed analysis of irreversible loss 
mechanisms. The effect of slow changes in the magnetic field upon the 
omnidirectional intensity of particles was considered and calculations 
for the ARGUS time period show that the changes are not inconsequential. 
Ry making conservative assumptions concerning the change in the magnetic 
field during the main phase decrease following Pass 414, it was shown 
that the relative intensity could change by 15% due to the betatron 
acceleration alone. 
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Table Cl (Ref. C2) Characteristics of G-M counters, Explorer 4 


Detector 


Shielding 


Sensitive 

to: 


Geometric 

Factor: 


D (Channel 1) 

"Shielded" 
Geiger-Mueller 
counter (Anton 
302) Cylinder 
approximately 
7 mm x 9 mm 
Scalar: 64 


1.2 gVcm Fe 
+1.6 gm/cm 2 
Pb (minimum) 


Electrons of 
E e > 5 MeVj 
Protons of 
E^40 MeV, 
X-rays of 
E> 80 KeV 
with low 
efficiency 


Omnidirectional 
geometric factor G : 
0.14 cm 2 for 0 
minimum stopping 
power, 0.82 cm 2 
for penetrability 
greater than 7 
gm/cm 2 


C (Channel 3) 

"Unshielded" 
Geiger-Mueller 
counter (Anton 
302) Cylinder 
approximately 
7 mm x 9 nun 
Scalar: 2048 


2 

1.2 gm/cm Fe 
(minimum) 


Electrons of 
\> 3 MeVj 
Protons of 
E l> 30 MeV, 
X-rays of 

E> 20 KeV 
with low 

efficiency 


Omnidirectional 
geometric factor G : 
0.14 cm 2 for 0 
minimum stopping 
power, 0.70 cm 2 
for penetrability 
greater than 5 
gm/cm 2 


C14 



t 


Table C2. Shielding and Range-energy data for G-M Detector C, D 

Detector 

Shielding 

Equivalent Al. 

Energy 


(gVcm 2 ) 

2 

(gVcm ) 

(MeV) 

C 

1.2 Fe 

1.46 

2.95 

D 

1.2 Fe 

+ 1.6 Pb 

5.89 

| 

7.6 
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4 
















Table C3. Explorer 4 Passes "over' - Huntsville, Alabama selected for 
omnidirectional radiation analysis 


ARGUS Event 1 (time of burst: 0230 UT August 27, 1958) 


Date 

Pass 

Pass 

Pass 

August 27, 1958 

414 

415 

4 16 

28 

427 



29 • 

440 



30 

453 

454 


31 

466 



September 1, 1958 

479 



2 

492 



3 

505 



4 

518 




ARGUS Event 2 (time of burst: 0320 UT August 30, 1958) 


Date 

Pass 

Pass 

Pass 

August 30, 1958 

453 

454 

455 

31 

466 



September 1, 1958 

479 



2 

492 



3 

505 



4 

518 
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2 

Pig. Cl. Geometrical factor, G q (cm ), versus the energy (MeV) 
necessary for transmission by the various thicknesses 
of materials. Circled points are from Table C2, using 
Eq. (Cl), assuming a linear dependence of G q with E 
between plotted points and constant values for larger 
energies (as shown). 
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Pig. C2. Omnidirectional count-rates vs time for a typical Pass of 

Explorer 4 in the ARGUS Event 1 region. The solid lines 

are average count-rates drawn through points not shown. 

Curve A. Unshielded G-M counter (Detector C, Channel 3) 
which detects electrons of energy greater than 
3 MeV. The background has been drawn in by 
estimation. 

Curve B. Shielded G-M counter (Detector D, Channel 1 ) 

which detects electrons of energy greater than 
7.6 MeV. The background has been drawn in by 
estimation. 

Curve C. Ratio of A/B, unshielded to shielded G-M count- 
rates. The ratios are of count-rates above the 
background. 

Curve D. Difference A - B, unshielded minus shielded G-M 
count-rates above the background. This curve 
shows the omnidirectional count-rates due to 
electrons between 3 and 7*6 MeV, a roughly 
"monoenergetic" count-rate curve. 


C18 





Pig. C3. 


Patn of satellite Explorer 4 in B-L coordinates, 
regions of ARGUS Events 1 and 2 are shown also. 


The 


C19 



(GAUSS) 







Fig. C4, 


Planetary magnetic three-hour 


range K 


P 


indices during 


the time interval of ARGUS Events 1 and 2. 
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sudden commencement 











Fig. C$. Logarithmic plot of true count-rate of the unshielded 
G-M counter (Detector C, Channel 3) versus the elapsed 
time since detonation of ARGUS Event 1. The solid lines 
are least squares fit of Eq. (C3) for the data of the 
sequence of Passes 414, 427, • •• 5*8. n = 1.13 for 
both curves. Intersection of the two straight lines 
with the ordinate axis yield values of C Q , Eq. (C3). 

The peak intensities were centered at L = 1.72. 
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Pig. C6. Count-rate ratio, unshielded {Channel 3)/shielded 

(Channel 1) omnidirectional G-M counters versus 

elapsed time T e since ARGUS Event 1. Both G-H 

count-rates were corrected for background count- 

rates. The dashed line is the plot of C = C t 

0 

Eq. (C3). 
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L=1.72 





Pig. C7. The geometrical width at half-maximum of the ARGUS 
Event 1 region versus the elapsed time T g since the 
event. (Pass 453 was the first to contain ARGUS 
Event 2 data.) 
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Fig. C8. True count-rate of the unshielded (Detector G, 

Channel 3) G-M counter versus the elapsed time 

T since ARGUS Event 2 for selected L shells, 
e 
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Pig. C9. Logarithmic plot of true count-rate of the unshielded 

G-M counter (Detector C, Channel 3) vs the elapsed 

time T since ARGUS Event 2 for selected L shells, 
e 

Passes 505 and 518 were during a period of high magnetic 
activity. 


C25 



(counts/ sec) 






PASS 



TO* 



(hours) 


<o 

*0 

<* 


1 



cs 

Woo 


<* 

V t** 

w 10 

1 

1 

1 1 



Pig. CIO. Count-rate ratio, unshielded (Channel 3)/shielded 
(Channel 1) G-M counter vs elapsed time T q since 
ARGUS Event 2. L = 2.11 is the region of peak 
intensity of the ARGUS Event 2 region. The numbers 
above the plotted points are satellite Pass (over 
Huntsville, Alabama) numbers. 
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’ 0 ° 10 1 10 2 

T e (hours) 


Pig. Cll. Critical energies below which the magnetic moment 

is conserved in interactions with hydromagnetic 

waves for particles mirroring at geomagnetic 

latitudes greater than 30° assuming the hm wave 

is cut off sharply at ^max. Shown are Hnax = 

1, 5, and V, sec" 1 (Ref. C19). 

1 * 


C2? 



^critical (MeV> 



L (earth radii) 


Pig. C12. Regions in space where whistlers may be effective 
in precipitating electrons into the atmosphere. 
The plot is from Eq. (C4) with the heart of the. 
whistler range being 1 to 10 kc/s (Ref. C20). 
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Pig. C13. Number of beta rays per fission from U per energy 

? 

interval for various energies E in units of m Q c . 

The data is taken from Heller (Ref. C2<f). 
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Fig. Cl4. Observed omnidirectional, true initial count-rate C Q 

contours vs B and L for ARGUS Event 1 (Detector C, 

unshielded G-M counter, Channel 3). C is the true 

o 

-1 13 

count-rate corrected for decay assuming the t 
law, Eq. (C3). The Initial time reference for the 
decay law correction is the beginning of ARGUS shell 
counting for Pass 427. The solid line count-rate 
contour curves are best fits to the data from Passes 
415, 4l6, 427, 440, 453, 466, 479 (points not shown). 
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Pig. C15. Omnidirectional flux density J Q contours vs B 

and L for ARGUS Event 1 (Detector C, unshielded 
G-H counter, Channel 3) obtained from the observed 
omnidirectional true count-rate C Q data (Pig. Cl4). 
The conversion is given by J q = 67 .C q . 




OMNIDIRECTIONAL FLUX DENSITY 



1.68 1.70 


D. Summary of M.S. Thesis, Joseph R. Ockersz, S.J., "Omnidirectional 
Flux Densities of Geomagnetically Trapped Radiation Explorer 4 
Pre-ARGUS Times"; Saint Louis University, Physics Department 
(May 1966). 

Dl. Introduction 

The purpose of the present Summary is to present experimental 
aspects of the trapped radiation, including the omnidirectional charged 
particle flux densities in the geomagnetic field, for pre-ARGUS event 
times of Explorer 4 Satellite data telemetered over Huntsville, Alabama. 
Such an investigation is important toward understanding the "background" 
radiation (natural radiation) against which the artificially injected 
ARGUS radiation is observed. 

D2. Experimental Details 

Instrumentation and Calibration of Explorer 4 
The University of Iowa group designed the detection instruments. 
Van Allen has published a complete description of the detection instru¬ 
ments (Ref. D2), Reference D3, Fig. 2,1, p. 16, shows a block diagram 
of the four detectors A, B, C and D, and the corresponding telemetry 
channels 2 and 5, 3 and 1. (See also Preface P2, Fig. P2-1). 

For the present thesis the omnidirectional count-rates were 
analysed, involving data obtained by G-M Detectors C and D. Preface 
P2 of this Report summarizes the characteristics of Detectors C and D. 
Table Dl lists the telemetry data available for the present analysis 
and Fig. Dl shows the location of Satellite Explorer 4 for some of the 
Passes. 

The count-rates of the G-M counters must be corrected for 
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"deadtimes". For detectors C and D, Anton 502 G-M co-’nters, the cqua 
tion for the correction is 


C = Ce^ C 
obs 


(Dl) 


where C , is the observed count-rate and C is the true or corrected 
obs 

count-rate and 7« , the deadtime, is 62.5_+1.3 microseconds. 

Further instrumentation and analysis details are found in DASA 
Report WT 1613, Ref. (D3). 


D3. Results 

A complete set of Figures showing the true count-rates vs time 
for the different B-L regions of the pre-ARGUS Explorer 4 times is 
found in Ref. Di. 

These Ref. Dl Figures are here compactly summarized by three new 
Figs. D2 (constant omnidirectional count-rate contours in B-L coordi¬ 
nates) , D3 (omnidirectional flux densities vs B), D4 (three dimensional 
plots of omnidirectional true count-rates vs B and L). 

The slope of Log(«jT) vs Log(B) straight line (derived from 
Fig. D4) compares closely with earlier published results (Ref. D4). 

In general it has been found (Ref. Dl) that Passes 2, 13, 26, 52, 78, 
107, 173, 207 and 367 recorded small flux densities, between 10 and 30 
particles cm'^sec"*; while passes 4, 43, 121, 272, 285, 298, 311, 324, 

337, 375 and 401 recorded large flux densities between 100 and 20,000 
-2 -1 

particles cm sec . Passes of small measured flux densities correspond 
to low altitude (250 < altitude < 700 km). Tne low count-rates are 
due chiefly to non-trapped cosmic radiation in space. The Passes of 
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larger measured altitudes (1000 altitude < 1500) are specified by 

0.17 < B < C.30, i.35 < L <1.65 and 0.35 < B <0.45, 3.5 < L ^ 5.5. 

The larger omnidirectional flux densities were recorded in the Van 

Allen region. 
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Fig. D1 Map shov/ing location of Satellite 1958€(Explorer 4) for 

Passes 2, 4, 13, 26, 43, 52, 63, 78, 107, 121, 172, 207, 272, 
285, 298, 311, 324, 337, 367, 375, 401 and 414 tracked by the 
U. S. Army Tracking Station at Huntsville, Alabama. Except 
for Pass 414, these Passes are for pre-ARGUS times. The 
Turber in parenthesis is the Pass number, the other numbers 
are the satellite altitudes at the beginning and end of the 
tracking period for the Pass. 





Fig. D2 Omnidirectional true constant count-rate contours in B-L 
coordinates. Data is for pre-A'IGUS times, shielded G-M 
detector D (channel 1) for the region defined by 
1.27 < L < 1.30 and 0.20 < B < 0.23. To obtain Jq, the true 
omnidirectional count-rates are multiplied by 8.4*, a reduc¬ 
tion constant derived from the calibration of the detector, 
(♦provisional, subject to correction) 
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Fig. D3 Omnidirectional flux densities vs magnetic flux density B. 

Data is for pre-ARGDS times, shielded G-M detector D (channel 1) 
for the region defined by 0.20 <• B < 0.23 and 1.27 < L < 1.30. 
Log J Q vs B is plotted instead of log J Q vs log B in order to 
obtain somewhat greater spread in data. The data may be repre¬ 
sented by the equation = k B n , whence n = 30, approximately. 
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Fig. D4 Three dimensional plot of omnidirectional true count-rates C 
(shielded G-M counter. Detector D, channel 1) (number sec"*) 
vs B and L. Sheets of constant L are shown intersecting the 
contours of constant count-rates. Data is the same as for 
Figs. D2 and D3. 





















Table Dl. Explorer 4 telemetry, pre-ARGUS times data available at 
Saint Louis University, Physics Department; Fasses are 
over Huntsville, Alabama Tracking Station. 


Pass 

No, 

Date 

Time 

1 

July - 26 

1658-1702 

2 

26 

1855-1858 

3 

26 

205^2057 

4 

26 

2240-2257 

5 

July - 27 

0036-0042 ! 

13 

27 

1650-1655 

15 

27 

1348-1851 

16 

27 

2042-2052 

17 

27 

2240-2247 

26 

July - 28 

1449-1451 

27 

28 

1642-1648 

28 

28 

1840-1844 

42 

July - 29 

2026-2035 

43 

29 

2220-2227 

52 

July - 30 

1430-1437 

53 

30 

1627-1632 

64 

July - 31 

1421-1426 

66 

31 

i6l8-l622 

67 

31 

1816-1818 

69 

31 

2205- 

70 

Aug. - 1 

0000-0007 

77 

1 

1220-1222 

78 

1 

1412-1420 

79 

1 

1609-1614 

80 

1 

1806-1811 

81 

1 

2001-2008 

82 

1 

2154-2208 

94 

Aug. - 2 

1951-1957 

104 

Aug. - 3 

-0841 

107 

3 

1941-1949 

108 

3 

2130-2141 

109 

3 

2331-2335 

121 

Aug. - 4 

2125-2134 

122 

4 

2321-232? 

133 

Aug. - 5 

1921-192? 

134 

5 

2113-2126 

147 

Aug, - 6 

2200-2210 

148 

6 

1257-1305 

169 

Aug, - ? 

1257-1305 

172 

Aug. - 8 

1845-1852 

173 

8 

2039-2050 

182 

Aug. - 9 

1243-1249 

189 

9 

1820-1828 


Pass 

No. 

Date 

Time 

199 

Aug. - 10 

2015-2025 

207 

11 

1025-1035 

208 

11 

1220-1228 

220 

Aug. - 12 

1010-1021 

221 

12 

1807-1815 

222 

12 

1405- 

22 + 

12 

1754-1801 

225 

12 

1948-1957 

233 

Aug. - 13 

0957-1008 

234 

13 

1153-1202 

237 

13 

1741-1748 

238 

13 

1935-1943 

246 

Aug. - 14 

0950-0953 

24? 

14 

1139-1149 ! 

250 

14 

1725-1734 

251 

14 

1923-1929 

259 

Aug. - 15 

0924-0938 

260 

15 

1125-1132 

264 

15 

1907-1915 

2? 2 

Aug. - 16 

0910-0920 

273 

16 

1107-1118 

2?6 

16 

1657-1705 

277 

16 

I852-I858 

285 

Aug. - 17 

0855-0906 

286 

17 

1056-1102 

298 

Aug. - 18 

0840-0854 

302 

18 

1627-1633 

310 

Aug. - 19 

0630-0642 

311 

19 

0823-0833 

312 

19 

1022-1032 

313 

19 

I6ll-l6l7 

324 

Aug. - 20 

0836-0809 

325 

20 

1004-1015 i 

328 

20 

1555-1602 

337 

Aug. - 21 

0?48-0806 

338 

21 

0949-0959 

341 

21 

1538-1544 

355 

Aug. - 22 

1521-1527 

367 

Aug. - 23 

1503-1510 

375 

Aug. - 24 

0510-0518 

401 

Aug. - 26 

0433-0440 

406 

26 

1408-1414 

413 

26 

2225-2259 
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E. Summary of Pertinent Analytic Sections of Ph.D. Dissertation, 

Joseph M. Palkeday, "Interpretation of Directional Flux Densities 
in ARGUS Shells, Explorer 4 Satellite Data"; Saint Louis Univer¬ 
sity, Physics Department (March 1966). 

El. Introduction 

The method used for obtaining the directional flux density J(0 1 ) 
-2 -1 -1 

(particles cm sec ster ; 0 ! , angle between scintillation counter 

axis P and J) follows that of Lundquist, Nauraann, and Weber (Ref. El) 

in the reduction of the integral equation relating the corrected direc- 

-2 -1 

tional count-rate C(0) (number cm sec ) and J(0') to matrix form 

A 

(0, angle between P and a plane perpendicular to the magnetic flux 
density vector B). The matrix elements S^ (Ref. E2, p. 13) associated 
with the response function have been evaluated to a greater accuracy 
resulting from a more precise definition in the present analysis of the 
limits of integration in the integral equation relating C(0) and J(0‘) 
(Ref. E2, Eq. (1.6)). It has been found that the kernel of this inte¬ 
gral equation has a vanishingly small, determinant and as a consequence 
established methods of solving the integral equation for J(0') values 
fail to produce physically meaningful results. The kernel of this 
integral equation is highly sensitive to experimental errors inherent 
in the C(0) data. The present method replaces the inverse of the 
kernel matrix by an approximate matrix (Ref. E2, Eq. (1.26)) which is 
insensitive to small errors in the counting rate data. 


El 


E2. Summary of the Analytic Procedure 

The complete details of the analytic procedure outlined in El 
above are given in the Paikeday Ph.D. Dissertation (Ref. E2/. The 
procedure involves the following steps: 

(1) Reduction of the integral equation to a matrix form (see Ref. E2 
PP. 7-13). 

(2) Representation of the directional scintillation detector (Detector 
A, Channel 2) response function S(a) (Ref. E2, Pig. 1.2, p. 149) 
in analytic form (see Ref. E2, Eq. (1.10)) to yield J(0*) values. 

(3) Inversion of the matrix equation (see Ref. E2, pp. 14-33). 

E3. Conclusions 

With the data available as discrete points of the C(Q) function, 
the solutions obtained by using the inverse of the S-matrix are found 
to be non-physical. The process of getting correct physical solutions 
does not in any way remove the ill-conditioning of the original matrix 
S representing the kernel of the integral equation but rather supplies 
information about the type of errors that enter the counting rate data. 

The method developed in the present analysis is found to yield 
meaningful J( 0 ') values for the available C( 0 ) data. 

E4. References 

El. Lundquist, C., Naumann, R., and Weber, A. H., "Directional Flux 
Densities and Mirror Point Distributions of Trapped Particles from 
Satellite 1958 Epsilon Measurements", J. Geophys. Res. 67 , 4125- 
4133 ( 1962 ). 
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